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EXPERIMENTAL AND ANALYTICAL INVESTIGATION OF 
THE FLUTTER OF FLAT BUILT-UP PANELS 
UNDER STREAMWISE INPLANE LOAD 

By Robert W. Hess 
Langley Research Center 

SUMMARY 

This paper presents both experimental and analytical flutter results of an investi- 
gation of the effects of streamwise inplane load on 13 panels at Mach numbers of 1.57 
and 1.96. The panels were supported by longitudinal stiffeners and transverse frames, 
the spacing of these members allowing a variation in the length -width ratio from 1.0 
to 4.2. The experimental results are compared with calculated results for six of the 
panels. Two sets of calculations were made by using linearized quasi-steady aerody- 
namic theory in conjunction with a Galerkin procedure. The calculations differed in the 
treatment of structural damping. One set of calculations employed calculated vibration 
mode shapes for an unstressed panel and the structural damping was applied to panel 
bending only. The second calculation employed calculated vibration mode shapes for the 
panel in a stressed condition, and structural damping was applied to both panel bending 
and compression. 


INTRODUCTION 

There are many variables that affect panel flutter and a major problem involved in 
the study of panel flutter is the isolation and determination of the significance of these 
variables. The flutter boundary for a particular panel is a complex function of such var- 
iables as panel configuration; midplane stress; edge restraint; flow angularity; dynamic 
pressure; Mach number; pressure differential and gradient; structural, aerodynamic, 
and acoustic damping; and boundary -layer thickness. The destabilizing effect of mid- 
plane stress has been the subject of many experimental and analytical investigations 
since the work of Sylvester and Baker (ref. 1). Reviews that contain some results of 
these efforts may be found in references 2 to 5. The agreement between experiment and 
theory has often been poor because of the inability to define or control adequately the 
variables affecting the onset of flutter. In addition, as pointed out in reference 6, the 
misapplication of damping in analytical solutions appears to have been the cause of much 
of the lack of agreement for stressed panels. 



The purpose of the investigation reported in the present paper was to examine the 
flutter of flat isotropic panels subject to inplane compressive loads in the streamwise 
direction only and constructed to simulate practical edge conditions. A series of 
13 multibay panels was constructed so that for each, the test panel was the center bay of 
a built-up structure. The center bay was supported by longitudinal stiffeners and trans- 
verse frames whose spacing defined test-panel length -width ratios ranging from 1.0 
to 4.2. The dynamic pressure and frequencies at flutter were calculated by two methods 
for six of the panels, and these analytical values were compared with the experimental 
results. Linearized aerodynamic theory for two-dimensional flow in which the aero- 
dynamic damping term was retained (also known as ’’modified piston theory”) was used 
in both methods in conjunction with a Galerkin (modal) procedure. Mode shapes calcu- 
lated with zero inplane stress were used in the first method in which the structural 
damping appeared only with the bending terms. Mode shapes calculated with inplane 
stress were used in the second method in which the structural damping was applied to 
both the bending and the inplane load terms. 

Appendix A presents the measured flow conditions over the panel. Appendix B 
discusses the techniques and apparatus used in measuring panel mode shapes, frequen- 
cies, and damping. Appendix C gives the details of the two analyses used in the present 
paper. 


SYMBOLS 

This list contains only those symbols that are significant to the entire paper. 
Symbols are defined in appendix C that are restricted in meaning to that section of the 
paper. 

Am>An complex modal amplitude coefficients 

aj speed of sound over panel 

B parameter defined by equation (Cll) 

Cq,Cj,C 2 modal integrals defined by equations (4) 

^mn^mn^mn modal integrals defined by equations (8) 

Eh D 3 

D panel flexual rigidity, — — 

12(1 - */ 2 ) 


2 



E 

f 

g 

h 

i = \Ti 


Young’s modulus 

panel frequency 

structural damping coefficient 

aerodynamic damping coefficient, equations (13) 

thickness 


k 

l 

M, 


M 


t,e 


N, 


- N x 
N = 
x ,2 
yl 


n t 


- N v c i 

N y — \ 


reduced frequency, ico/V^ 
length of test panel 
local Mach number 
torsional stiffness of stiffeners 

inplane loading in stream direction, positive in compression 


inplane loading in cross-stream direction, positive in compression 


Ap static pressure differential 

q dynamic pressure 

R radius 

t time 

v i 


velocity of airstream over panel 



w 


width of test panel 


x,y,z 


orthogonal coordinates (see fig. 11) 


streamwise asymmetric deflection function 


X c 


streamwise symmetric deflection function 


spanwise deflection function 


panel deflection 




eigenfunction in arguments X^,Xg 


y 


mass per unit area 




nondimensional coordinate in y-direction, 2- 

J ’ w 


dynamic -pressure parameter, 

j3 M D 

density ratio. -2— 

P x l 

Poisson's ratio for isotropic plate 


b 

pi 


nondimensional coordinate in x-direction, - 
air density over panel 


inplane stress in stream direction, positive for compression 


cr 


Q 


static buckling stress with no airflow 


flutter eigenvalue parameter, (co r / co)^ ( 1 + ig) 


co 


circular flutter frequency 
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circular reference frequency, 


w r 

Subscripts: 



circular natural frequency mode m 


m,n mode numbers 

p panel 

ad aerodynamic damping 

r reference 

s stiffener 


APPARATUS 
Wind Tunnel 

The flutter tests were conducted in the Langley Unitary Plan wind tunnel at free- 
stream Mach numbers of 2 and 1.6. The dynamic pressure in that tunnel is continuously 
variable with maximum levels of 2550 psf (122 kN/m 2 ) and 2400 psf (114.9 kN/m 2 ) at 
free-stream Mach numbers of 2 and 1.6, respectively. 

The splitter plate on which the panels were mounted (fig. 1) projected into the air- 
stream 15 inches (38 cm) from the tunnel sidewall so that the flutter tests could be con- 
ducted free of the tunnel boundary layer. To maintain the smallest possible pressure 
gradient on the surface of the panel, the splitter plate had a negative angle of attack 
of 1.05° with respect to the free stream. The resulting local Mach numbers over the 
panel were 1.96 and 1.57 for the respective free-stream Mach numbers of 2.0 and 1.6. 
The local static and dynamic flow conditions over the panel are given in appendix A. 

Vacuum Chamber 

The panel vibration characteristics were determined in the 55 -foot-diameter 
(16.764 m) vacuum chamber of the Langley dynamics research laboratory. The vacuum 
system is capable of evacuating the chamber from atmospheric pressure to 0.2 mm 
mercury absolute. 
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Splitter Plate and Loading Mechanism 

The splitter plate (fig. 2), which had a cavity 1.8 inches (4.572 cm) deep, 

38.42 inches (97.587 cm) long, and 24.38 inches (61.925 cm) wide, contained the 
hydraulic mechanism for applying inplane loads to the panel. The leading edge of the 
panel was bolted to the splitter plate and the trailing edge was fastened to a movable 
mounting plate. A guide plate was supported from the back of the splitter plate cavity 
by three low-friction guides. The guide plate constrained the mounting plate to recti- 
linear movement only. Inplane compression loads (8000 lb (35 586 N) maximum) were 
applied at the trailing edge by means of hydraulic pressure acting on eight pistons which, 
in turn, loaded the movable mounting plate at the trailing edge. A single piston was used 
to return the system and to apply tension to the panel. 

The streamwise edges of the built-up panel could be clamped to the splitter plate 
(fig. 1) by clamps which, upon being released, relieved accumulated shear stress due to 
inplane loading. Each clamp was pulled down by the vertical movement of five recessed 
bolts. The movement of these bolts (clamped or released) was actuated by remotely 
controlled pneumatic jacks. The leading edge of the clamp was bolted to the splitter 
plate and a flex hinge, which had been milled into the clamp ahead of the panel leading 
edge, allowed translation and rotation of the remaining portion of the clamp. 

All openings to the splitter plate were sealed with appropriate O-rings, gaskets, 
or sheet rubber. Pressure in the cavity was controlled by bleeding in or evacuating air 
through a 1-inch (2.54-cm) tube. 


INSTRUMENTATION 

The basic instrumentation consisted of four four-arm strain-gage bridges mounted 
on the test panel. Three of these measured compression; one at the leading edge, one at 
the trailing edge, and one next to the middle of a longitudinal stiffener at one side of the 
test panel. The fourth strain-gage bridge measured bending at the trailing edge of the 
panel. Its dc output was monitored on a digital voltmeter, and its ac output was moni- 
tored on an oscilloscope. Two pressure cells (5 psi (34.475 kN/m2)), which measured 
the pressure differential between the cavity behind the panel and the surface of the panel, 
were used to establish a minimum pressure differential. 

The output of the compression gages, panel pressure cells, and the pressure cells 
used for measuring the hydraulic pressure of the loading apparatus was monitored on 
self-balancing potentiometers. These outputs and that of the bending gage were also 
recorded by a recording oscillograph and a frequency -modulated tape recorder. 
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Thermocouples were mounted at the center of the panel and at the leading edge 
adjacent to the transverse frame. The temperatures at these two locations were read on 
recording potentiometers. 


PANELS 

Configuration 

The panels were built-up structures, as shown in figures 3 and 4, in which only the 
center bay constituted the test panel. The method of construction was chosen to provide 
a realistic edge restraint for the panels. The length and width of the center bay (mea- 
sured from rivet line to rivet line) was varied by using different spacing of the trans- 
verse frames and of the principal longitudinal stiffeners (those adjacent to the center 
bay). The principal longitudinal stiffeners were continuous and welded to the doublers 
at the leading and trailing edge of the structure whereas the auxiliary stiffeners were 
interrupted by the transverse frame. The auxiliary stiffeners were incorporated to 
raise the buckling loads and the flutter dynamic pressure of the adjacent panels to levels 
higher than that of the test panel. In addition, some bays adjacent to the test panel were 
stiffer since they were constructed with a thicker skin which overlapped the test panel at 
the stiffener. This construction was used only for high length -width ratio panels. A 
summary of the panel configurations is given in figure 5, and a table of the physical 
properties is given in table I. Two orifices were located in the skin of the built-up panel, 
one ahead of the forward transverse frame, the second just back of the aft transverse 
frame, to provide a means to measure the pressure differential between the cavity and 
the surface of the panel. 


Buckling Stress 

The ratio of inplane stress to static buckling stress o/a cr is a primary variable 
in this investigation. The buckling stress was calculated and determined experimentally 
for each panel. The calculated values were determined from the equations of refer- 
ence 7. The torsional stiffness M t q of the edge supports was calculated by using 
shear flow analysis. (See ref. 8, for example.) In order to verify the accuracy of the 
calculated values of Mj. q , seven test specimens such as shown in figure 6 were con- 
structed and tested, and the important dimensions are given in table II. Experimental 
values of Mj- q were between 20 and 27.7 percent lower than the calculated values. 

On five of the specimens, the rivet spacing was reduced by a factor of two and, as a 
result, the difference between calculation and experiment was reduced to between 2.7 
and 15 percent. The design rivet spacing was not changed on the test panels, but the 
calculated values for the torsional stiffness were consequently reduced by 24 percent as 
given in table I under "Adjusted torsional stiffness." 
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Experimental values of buckling stress are, in general, dependent on the means of 
loading (by displacement or by force), panel support, and initial deformation. The neces- 
sity of using pneumatic clamps to secure the side edges of the built-up panel to the 
splitter plate imposed an additional restraint and all buckling measurements were made 
with the clamps in the clamped condition. The load was applied in increments and the 
clamps were cycled (released and clamped) at each increment. In all cases the actua- 
tion of the clamps increased the inplane stress somewhat. This condition was undoubtedly 
due to the fact that the edges of the panel did not precisely conform to the splitter plate 
along the length of the panel edge and the final increment of stress was actually applied 
by the side clamps. 

Zero stress was consistently obtained by initially putting the panel in tension, 
releasing the hydraulic pressure, and then repeatedly cycling the clamps until no change 
in strain was measured in the clamped condition. This procedure gave repeatable values 
of buckling stress and repeatable values of zero strain after compression loading. 

Table I lists two values of experimental buckling stress because in most cases, the mea- 
sured stress at buckling differed between the leading-edge (L.E.) and the trailing-edge 
(T.E.) strain -gage locations. In most cases, the location with the highest stress level 
had a higher stress -loading slope from no load to. buckling as is shown in sketch (a). 

The thickness of critical members — skin, stiffener, and frames — differed on 
panels 7 to 13. The requirement that the stiffeners be stable beyond the critical buckling 
stress of the test panel resulted in large values of stiffness ratio gy/D and essen- 
tially clamped-edge conditions. 


Frequencies 

Panel frequencies were measured at atmospheric pressure before and after the 
flutter tests. The small cavity depth of 1.8 inches (4.57 cm) in the splitter plate pre- 
cluded meaningful comparison with calculated frequencies because of the stiffening and 
damping effects of the air in a shallow cavity. These effects were explored experiment- 
ally and analytically in references 9 and 10, respectively. At the completion of the flutter 
tests, the splitter plate was installed in the vacuum chamber and the natural frequencies, 
mode shapes, and structural damping were measured for panels 1 to 5 at levels of inplane 
load and ambient pressure that corresponded to those of the experimental flutter bound- 
ary. The details of the apparatus and procedure are described in appendix B. 

For comparison with experiment and for use in the flutter analysis, the natural 
frequencies and mode shapes were calculated for panels 1 to 6 by using the equations of 
reference 7. Figure 7 gives the calculated natural frequencies as a function of inplane 
stress. Figure 7 also contains the corresponding measured frequencies for panels 1 
to 5. For comparison purposes, the experimental values at zero stress are given at the 
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lowest density for which they were obtained. A comparison of measured and calculated 
mode shapes is shown in figure 8 for the first and third modes of panel 3 at stress 
ratios o/a cr of 0, 0.20, and 0.45. These mode shapes were measured after flutter 
testing and are considered to be typical. The greatest difference between calculated and 
measured mode shapes was in the lower modes, and this difference was amplified by 
increased inplane load. It should be noted that all symmetric modes were excited by a 
shaker at the center of the panel. 
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The solution of the equations of reference 7 for the panel frequencies and mode 
shapes as a function of inplane stress was programed by using Newton iteration proce- 
dures. Equations (16) to (19) of reference 7 are of two types depending on panel geom- 
etry and inplane load. For zero inplane stress, the eigenfunction is composed of trigo- 
nometric and hyperbolic functions. As the compression load is increased, the solution 
enters a region (shown as dashed lines in fig. 7) in which the lower eigenfunctions are 
composed only of trigonometric functions. As is found in reference 7, there are levels 
of stress in this region where with changing stress, panel modes evolve so that inflec- 
tion points and internal node lines appear and disappear. During a preliminary test of 
the mode-shape surveying apparatus in the vacuum chamber, before the flutter tests, 
modes such as shown in figure 5 of reference 7 were encountered. However, because 
of the difficulty in mapping the mode shapes of a damped system at high stress levels, it 
was decided to restrict measurements to conditions of inplane loading below a stress- 
ratio cr/a cr of 0.55. 

The frequencies measured at atmospheric pressure and with zero inplane load 
were lower after the flutter tests than they were before. With the exception of panel 3, 
the measured values of the panel frequency in a vacuum are significantly less (10 to 
20 percent) than the calculated values. This result coupled with the fact that the mea- 
sured buckling stress prior to flutter testing was close to the calculated value (except 
for panel 4) is a measure of the loss of stiffness caused by the flutter testing. 

Damping 

For no airflow, the damping of the panel vibration was affected by the interdepen- 
dent variables of cavity size, ambient pressure, inplane stress, mode, frequency, ampli- 
tude, and energy dissipation at structural joints. The ambient pressure and stress ratio 
were varied during the damping measurements of the various modes. Because of the 
number of variables and their interdependence, a plot of damping as a function of fre- 
quency results in a wide scatter band and no general conclusions are justified. The test 
results are tabulated in table III and were obtained in the vacuum chamber for the lower 
modes of panels 1 to 5. The damping coefficient g was determined from measurements 
of the vibration decrement. Most of the data fall between values of g of 0.015 and 0.05 
which is in the range cited in the published literature; that is, 0.01 for simple structures 
to 0.08 for complex, built-up structures. (See ref. 11.) 

It was anticipated that panel 5 might have higher values of damping because of its 
construction; that is, the center bay had a lap joint with the thicker adjacent panels. As 
can be noted in table III, the level of damping does appear to be higher than that of 
panels 1 to 4. It was not anticipated that an application of an inplane load would result 
in a decrease in measured damping, but this decrease did occur for modes 3, 4, and 5 
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of panel 5. A close look at the damping of panels 1 to 4 indicates isolated examples of 
this behavior. A possible cause for this behavior, which was not explored, may be the 
coupling of modes between the test panel and the adjacent bays. As was discussed by 
Ungar and Scharton (ref. 12), indirect excitation of one panel or system by another 
causes an energy exchange and higher recorded values of damping. The degree of 
coupling between bays would be expected to vary with inplane stress since the rate of 
change of the natural frequencies with changing inplane load would not be the same for 
the center bay and the adjacent bays. 

FLUTTER TESTS 
Procedure 

The general procedure was to set the tunnel flow to the desired dynamic pressure 
while maintaining the pressure in the cavity at a minimum difference from the pressure 
on the surface of the panel. The output of the two thermocouples was then monitored 
until both gave the same constant reading. The procedure for determining an effective 
zero-pressure differential was as follows: The inplane compressive load was incre- 
mentally increased until flutter was observed, at which time the differential pressure was 
slowly varied until maximum flutter amplitude was observed. The readings of differen- 
tial pressure from the two transducers at maximum flutter amplitude were used as the 
new value of zero differential pressure. 

After an effective zero differential pressure was established, the inplane load was 
released and the zero reference stress was recorded. Inplane load was then slowly 
increased incrementally. At each increment the clamps were cycled (released and 
clamped) to relieve any accumulated shear stress at the side edges. When the flutter 
boundary was approached in this manner, the first indication of flutter was the appear- 
ance of a small harmonic signal, superimposed on a response signal generated by aero- 
dynamic noise, which usually increased slowly in amplitude with time. In some cases it 
was necessary to increase the inplane load beyond the initial point before the amplitude 
increased. 

This procedure was not used for the case in which the minimum -dynamic -pressure 
flutter boundary was determined for a panel stiffened by the pressure differential. It was 
not prudent at the higher levels of dynamic pressure and differential pressure to release 
the clamps. To do so would result in the loss of stabilizing pressure differential and the 
possible loss of the strain gages or the panel because of violent flutter. 

At low dynamic pressures, pressure differential and inplane load could be applied 
and the clamp cycled without risk to the panel. The dynamic pressure was then raised 
until flutter was observed. At this point both dynamic pressure and inplane load were 
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slowly varied until the minimum dynamic pressure at flutter was obtained. The dynamic 
pressure was then dropped and the procedure repeated at a higher pressure differential. 

Experimental Results 

The experimental results at Mach numbers of 1.96 and 1.57 are presented in fig- 
ure 9 for flutter frequency and flutter dynamic pressure as a function of stress ratio. 

The stress ratio cr/a cr used in figure 9 is an average value. As was discussed in the 
section dealing with buckling stress, the buckling stress was not the same at each bridge 
location. At a given inplane load, the stress ratio was determined by using the buckling 
stress (J cr for each panel compression bridge location and taking the average of the 
resulting ratios. 

The data of figure 9 are characterized by a rapid decrease in the flutter dynamic 
pressure with increasing inplane stress. With the exception of two panels, the minimum 
dynamic pressure occurred at inplane stress ratios as low as 0.6. This result is in con- 
trast to the flutter behavior of thermally stressed panels which have a nearly equal load 
in the streamwise and cross-stream direction. For this type of loading (ref. 6), the 
minimum flutter dynamic pressure is at the transition point (the intersection of the 
buckled and unbuckled panel boundary). 

The flutter boundaries fell into two general categories at both Mach numbers: 
those which form a single smooth curve with increasing inplane stress, and those which 
exhibit abrupt changes. These changes are also accompanied by a jump in the panel 
flutter frequencies. An abrupt change in frequency can also be noted at or near zero 
stress for some panels whose boundaries fall in the first category. This flutter was at a 
frequency higher than that of most of the data, and the amplitude at a given stress ratio 
increased more slowly with time. Panel 12 (fig. 9(1)) is an example of the most extreme 
change in flutter boundary in that the high-frequency flutter was completely separate 
from the main body of the data and could be approached from a high and low level of 
stress ratio. The occurrence of high-frequency flutter in most cases prevented a good 
determination of the low-frequency flutter boundary near the zero stress level. 

The M^ = 1.57 data are at the same or higher dynamic pressure, at a given 
stress ratio, than the M^ = 1.96 data. This result is contrary to the expected result 
that the flutter dynamic pressure decrease with decreasing Mach number. This result 
may be accounted for by considering the effects of local flow conditions at the two Mach 
numbers. The aerodynamic noise and the static -pressure variation along the center of 
the panel were both greater at M^ = 1.57 than at Mj = 1.96. The panel response to 
the aerodynamic noise had the effect that the onset of flutter was more difficult to detect 
at a Mach number of 1.57 and the greater static -pressure variation had a stiffening effect 
on the panel. 
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The stabilizing effects of maintaining a positive pressure differential (internal 
pressure greater than static external pressure) on an inplane -loaded panel were inves- 
tigated with panel 11 at = 1.96. The results are given in figure 10 for flutter fre- 
quency and dynamic pressure as a function of static -pressure differential. Each point 
represents the minimum flutter dynamic pressure for the given differential pressure. 
Because of the stiffening effects of pressure differential, the inplane load at the minimum 
points increases with increasing Ap. The variation of the minimum dynamic pressure 
at flutter with pressure differential is relatively small near zero, but increased rapidly 
above a differential pressure of about 0.2 psi (1.379 kN/m2). At a Ap of 0.43 psi 
(2.96 kN/m^), the dynamic pressure at flutter is increased by a factor of eight; this 
increase is equivalent to a factor of two increase in panel thickness. 

ANALYSIS 

Based on the results of previous investigators, it is evident that structural and 
aerodynamic damping should be included in an analytical panel flutter model. Neglecting 
damping completely, as shown in reference 13, results in a flutter boundary (q as a 
function of cr/a cr ) that goes to an unrealistic zero dynamic pressure at points of fre- 
quency coalescence for panels with length-width ratios greater than 1.0. The addition of 
aerodynamic damping improves the results somewhat; however, an exaggerated decrease 
in flutter dynamic pressure is still conspicuous at points of frequency coalescence. 

Structural damping (that is, damping in harmonic motion which is proportional to 
amplitude and in phase with velocity in contrast to a viscous representation), in addition 
to aerodynamic damping, is used in two analytical models in the present paper. In the 
first model the structural damping coefficient modifies only the bending term. Calcu- 
lated mode shapes and frequencies for an unstressed panel are used in the calculation 
which is hereafter called the zero-stress-mode solution. The advantages of this model 
are demonstrated in reference 6. In the second model, structural damping is assumed 
to be equally effective in both bending and compression as assumed in references 13, 14, 
and 15. Calculated mode shapes and frequencies for stressed panels are used in this 
calculation which is hereafter called the stressed -mode solution. 

Flutter boundaries were calculated for panels 1 to 6 for a Mach number of 1.96. 

The panel properties are given in table I and the coordinate system is shown in figure 11. 
The aerodynamic loading was represented by linearized, two-dimensional aerodynamics, 
the aerodynamic damping term being retained (also known as "modified piston theory"). 

A brief outline of the two analyses, in which a Galerkin procedure is used, is presented. 
The integrals resulting from the derivation are given in appendix C. 
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Flutter Solution With Zero -Stress Modes 

The differential equations of motion for small deflections of an isotropic panel sub- 
ject to inplane loads with one side exposed to a supersonic airstream may be written as: 


( a 4 z 2 a 4 z a 4 z \/i . * N x a 2 z N y a 2 z 

W 4 + j2 w 2 H 2; jq 2 W 4 S „4) (1 + lg) + t 2 g^2 + w 2 dj] 2 


+ 



Vi_/az + Vi8z\ 

Pl %v 8t 1 a iJ 


(i) 


where 

l = r 

g structural damping coefficient 

N x ,Ny inplane loading in the x- and y -directions, positive in compression 

% = \/ M i - 1 

y mass per unit area of panel 

p ^ air density 

velocity 

The mode shapes that satisfy the boundary conditions for clamped edges with zero 
translation and with no airflow were calculated by using the equations of reference 7. 

An approximate solution to equation (1) was obtained by a Galerkin procedure where 
the solution was taken to be of the form 

N 

Zfetj.t) * Y, A n x n(£) Y t (n) e iwt (2) 

n=l 
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where X n is the nth mode shape in the x -direction, Yj is the assumed mode shape 
(one half-wave) in the spanwise direction, and the coefficient A n is, in general, com- 
plex. Substitution of equation (2) into equation (1), multiplication by Yj, and integration 
of the product with respect to 77 over the panel width results in the following ordinary 
differential equation: 


f AnV/ D(±X^ + p 


C i^n . 1 c 2 


2 W 2 Cq 

+ 5 ^ Xn " “ 2yXn + Pl lit ( iMXn + T X ") 


X " + ^^ X fl + le) + TT X n 


N x ^n 
€ 

V 


= 0 


(3) 


where the Roman superscripts on X n denote derivatives with respect to its argu- 
ment £. By following the notation of reference 7, 



and the Roman superscripts on Y 1 indicate differentiation with respect to its 
argument 77 . 

The differential equation for the free vibration of a panel with no inplane load con- 
sistent with the assumptions and notation used previously can be written as: 




+ 


2 Cl V 11 
l 2 w 2 Co n 


+ 



n ^n^n 


(5) 


which expresses the relationship between the elastic restoring forces and the normal- 
mode mass-inertia forces. 

Substituting equation (5) into equation (3) is valid since the modes and frequencies 
of equation (5) are panel normal modes. The resulting equation is 


£ w / [a£(l + ig) - u*]yX n + ^ x n + 


N v C 1 V l/ V 1 I 

i^ Xn+Pi %r Xn+ ^ Xn ' 


An = 0 


( 6 ) 
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Continuation of the Galerkin procedure by multiplying equation (6) by X m (£) and 
integrating over the length results (when m is successively 1,2,3,...N) in a set of N 
simultaneous equations in matrix form: 



C 1 N v 
C 0 yw 2 w 2 



nc (4) 

^x'-'mn 


,( 5 ) 

'mn 


yZ 2 Ci> 2 




M 



- 0 


where 


/ n = 1,2,3,4,...N;\ 

\m = 1,2, 3, 4,. ..N / '- 1 ' 


k = — 





1 


.(3) _ 

'mn 

C 1/2 

U/2 X “ Xnd| 

(=0 if m ^ n) 



,(4) _ 

mn ~ 

pl/2 jj 

[ 1/2 d « 


> 

(8) 

(5) _ ' 

mn - 

fl/2 j 

J_ 1/2 X m x n d « 

> 




For a nontrivial solution the determinant of the coefficients A n must equal zero. 
As is discussed in appendix C, the determinant was first put into a lambda form since the 
frequency term 1 /a> 2 appears in off-diagonal elements and then solved by means of a 
complex eigenvalue routine. 


Flutter Solution With Stressed Modes 
The differential equation for the stressed mode solution is: 


DA 


1 a 4 z 


a 4 z 


af 


Z 2 w 2 d^drj 2 


a^z 

W 4 ar? 4 y 


N x a 2 z 
i 2 a^ 2 


N y 9 2 Z 


w 2 a?^ 


(1 + ig) + y 


a 2 z 

at 2 


= .p 1 lL(lz + Xliz'\ (9) 
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Equation (9) is similar to equation (1) except that in equation (9) the structural damping 
for inplane compressive vibration is assumed to be equal to that for flexural vibration, 
both being multiplied by 1 + ig. 

Equation (9) was reduced to an ordinary differential equation as described in the 
previous section. The differential equation for the free vibration of a panel subject to 
inplane loads can be written as: 


A n°fe 4 IV) t J2r X ?*ir X » t I 5! »' t ^ X ») = r "» A n X n (10) 

u* w^ l w^ L, 0 / 


The ordinary differential equation becomes 


£ W |j^(l + ig) - w 2 ]yX n + Pi ^ iwX n + A n = 0 


( 11 ) 


Multiplying by X m (|) and integrating along the length results in the following N 
simultaneous equations in matrix form that express the dynamic equilibrium of flutter: 




^l k % 


’(5) 

•mn 


( A n> 


= 0 


/ n = 1,2,3,...N;\ 
\m = 1,2,3,...N / 


( 12 ) 


and 


a = 



(1 + ig) 


(jO 


r 



DISCUSSION OF ANALYTICAL FLUTTER RESULTS 
AND COMPARISON WITH EXPERIMENT 


The flutter boundary calculations of this investigation were highly dependent on 
both structural and aerodynamic damping. Retention of the aerodynamic damping term 
had a stabilizing effect. The ratios of the flutter dynamic pressure calculated with aero- 
dynamic damping to the calculated result without damping, and the corresponding frequency 
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ratios are given as functions of length-width ratio in figure 12. Structural damping and 
stress ratio were both zero for the calculations which were made by use of the zero- 
stress mode equations. Since aerodynamic damping is frequency dependent and since 
the natural frequencies of the panels increase with increasing length -width ratio (fig. 7), 
an increase in aerodynamic damping with length-width ratio is not surprising. Some 
insight into the large variation of the results in figure 12 may be obtained by considering 
the aerodynamic damping coefficient g a utilized in reference 16 for a semi-infinite 
panel: 


_ p l a l to 

S a yCO <y r 


(13a) 


which can be written as 


ga = 


7 T^\Jl>y 


(13b) 


The ratio of g a for the panels to that for panel 1 for a constant dynamic pressure 
is given in figure 13 for the panels in figure 12. Since the lengths of the panels are the 
same with the exception of panel 1, the variation of g a among the other four panels in 
figure 13 is due to panel thickness in D and y. 

The calculated flutter dynamic pressure for zero aerodynamic and structural 
damping agreed with the values computed from reference 17. The convergence of the 
present solutions was tested by increasing the number of modes used in the solution 
until the calculations gave small differences in the dynamic pressure and flutter fre- 
quency with increased number of modes. Eight modes were used for panels 1, 2, and 3; 
10 modes for panel 4; and 16 modes were used for panels 5 and 6. The results of ref- 
erence 17 were determined from a closed-form solution to the panel flutter equations 
and are based on the assumption of a spanwise mode, no damping, and piston-theory 
aerodynamic loading. 

Three flutter boundaries were calculated for panels 1 to 6 and are compared with 
the experimental data for = 1.96 in figure 14 in terms of flutter dynamic pressure 
as a function of stress ratio. The calculated and experimental flutter frequencies are 
also given as a function of stress ratio in figure 14. Two of the flutter boundaries were 
determined by using a structural damping coefficient g of 0.05 in the stress-mode and 
zero-stress -mode calculations. A structural damping coefficient of 0.05 was used since 
it was the upper limit for the bulk of the damping data. The third boundary was deter- 
mined by using zero structural damping in a zero-stress-mode calculation. The aerody- 
namic damping term was retained in the calculations of all three boundaries. 

The flutter boundary for the case of zero structural damping is shown as a dashed 
line in figures 14(a) to 14(f). For panels 2 to 6, the flutter boundary is similar to a 
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boundary in which aerodynamic damping is also neglected in that the boundary goes 
nearly to zero at points of frequency coalescence. This tendency is not noted in the 
experimental data. 

The flutter boundaries for the two calculations in which structural damping was 
retained are the same at low values of stress ratio. For panels 2 to 6 structural 
damping is stabilizing. As the stress ratio is increased, the two boundaries separate. 
The boundary from the stressed-mode solution (long dash) becomes irregular for panels 
with length-width ratios greater than 1.0. The zero-stress-mode boundary (solid line) 
remains smooth for all values of stress ratio and except for panel 1, the dynamic pres- 
sure at flutter is greater beyond the point where the boundaries separated. At the 
higher values of stress ratio (approaching 1.0), structural damping becomes less effec- 
tive in the stressed-mode solution and for panels 2 to 6 the results are essentially the 
same as those with a structural damping of zero. In contrast to the flutter dynamic 
pressure, the flutter frequencies from the two sets of calculations were in good 
agreement. 

The calculated results for panel 1 (Z/w = 1) are substantially different from the 
calculated results of the other five panels, in that the addition of structural damping was 
destabilizing for both types of calculations. This result is not uncommon in lifting - 
surface flutter (ref. 18) or in panel flutter (ref. 19). The effect of structural damping 
was more pronounced in the zero-stress-mode calculations. The experimental data of 
panel 1 fall uniformly below the calculated values of flutter dynamic pressure over the 
whole range of stress ratio a/cr cr in contrast to the results of the other five panels. 

The experimental data of panels 2 to 6 have the same general trend in the variation 
of dynamic pressure with a/ cr cr as the calculated data and the minimum experimental 
levels are in fair agreement with the zero-stress-mode calculations for g = 0.05. How- 
ever, there are some large differences between the calculated and experimental results 
at low levels of a/cr cr where the experimental results are confused by the occurrence 
of high-frequency flutter. 

With the exception of the high-frequency flutter, the bulk of the measured flutter 
frequencies were in fair agreement with the calculated values up to a stress ratio of 
about 0.70. Beyond this stress ratio the calculated flutter frequencies decreased more 
rapidly than the experimental data and go to zero above stress ratios of 1.0. 

A comparison of experimental and calculated results (zero -stress -mode shape 
calculation with a structural damping coefficient of 0.05) for M-^ = 1.96 is given in 

figure 15. The results are given in terms of the parameter xV3 ( where 

x l/3 w _ / 2 gw 3 

l 


l/3> 


as a function of length -width ratio for two conditions. The 
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calculated boundary for the minimum dynamic pressure with inplane streamwise load is 
shown as a solid line. The calculated boundary for zero inplane load is shown as a 
dashed line. The experimental data symbols are coded by panel thickness. 

As was discussed in the analysis section, a damping coefficient of 0.05 was chosen 
as an upper bound since the bulk of the measured damping coefficients fell below that 
value. The calculated results for the minimum flutter dynamic pressure reflect this 
choice in damping coefficient in that they are an upper bound to the experimental flutter 
data. A better fit of the calculated values to the experimental data would have been 
obtained if measured damping coefficients at all stress ratios had been used. 

The high-frequency flutter confused the experimental results near a stress ratio 
of zero. The limited experimental data at <j/cr cr = 0 shown in figure 15 that are com- 
patible with the linear theory used in the present paper are the result of extrapolation. 
The region of the high-frequency flutter at cr/a cr = 0 is shown as a shaded area in 
figure 15. 

CONCLUSIONS 

An experimental and analytical study has been conducted of the flutter character- 
istics of flat panels that were constructed with practical edge conditions. Thirteen 
panels were subjected to streamwise inplane loads and the study covered Mach numbers 
of 1.96 and 1.57. The flutter boundaries for six of these panels were calculated by using 
two-dimensional quasi-steady aerodynamics. Aerodynamic and structural damping were 
included in the formulation of the problem which was solved by use of a Galerkin pro- 
cedure. Structural damping was considered in two ways: effective in bending only, and 
effective in both bending and inplane compression. The investigation leads to the fol- 
lowing conclusions: 

1. The dynamic pressure at the onset of flutter decreased rapidly by 80 to 90 per- 
cent with increasing streamwise inplane compressive stress. For panels with length- 
width ratios greater than 1, the minimum flutter dynamic pressure occurred at com- 
pressive stresses that were somewhat less than the buckling stress, and the dynamic 
pressure subsequently increased when the panel was loaded beyond the buckling stress. 

2. Damping, both aerodynamic and structural, was found to have a pronounced 
effect on the calculated flutter boundary. 

3. The agreement between analysis and experiment depended on structural damping 
and the manner in which it was introduced into the formulation of the problem. Calcula- 
tions in which structural damping was introduced only in the panel -bending terms 
resulted in good agreement with experimental trends. 
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4. Limited examination of the effects of static -pressure differential acting across 
the panel indicated that static -pressure difference increased the inplane load required 
to cause buckling and resulted in an increase of the minimum flutter dynamic pressure 
by up to a factor of eight for the particular panel tested. 

5. Natural frequencies, mode shapes, and structural damping were measured for 
five of the panels at levels of inplane stress and ambient pressure that corresponded to 
the flutter condition. It was found to be impractical to map with sufficient accuracy and 
detail the number of mode shapes required for flutter calculations over the range of 
inplane stress of the flutter tests. 

6. Measured structural damping varied over a range from about 0.012 to 0.05 for 
variations in panel geometry, mode, inplane stress, and ambient pressure. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 20, 1969. 
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LOCAL FLOW CONDITIONS 

Local flow conditions were measured with an 0.5 -inch -thick (1.27 cm) calibration 
plate which was machined at the edges to fit the cavity so that its surface was flush with 
the splitter plate. Eighteen static -pres sure orifices were located 2 inches (5.08 cm) 
apart and 3 inches (7.62 cm) below the stream-alined center line of the plate. A 1-psi 
(6895 -N/m^) pressure gage was used to measure fluctuating pressures and was located 
6 inches (15.24 cm) aft of the vertical center line and 3 inches (7.62 cm) above the 
stream-alined center line. Provision was made to mount a boundary -layer rake on the 
center line at the leading and trailing edge of the plate. 

The pressure distributions for two levels of dynamic pressure at local Mach num- 
bers Mi of 1.96 and 1.57 are shown in figure 16. No values of Ap are given at the 
trailing edge of the plate because this area was in the wake of the boundary -layer rake. 
Pressure spectra, from 0 to 600 Hz, and overall levels expressed as ratios to dynamic 
pressure squared p rms are shown in figure 17. 

The pressure differential measurements in figure 16 show that at a local Mach 
number of 1.57, pressure differentials are two to three times greater than at a local 
Mach number of 1.96. The fluctuating pressure levels (fig. 17) were higher at a Mach 
number of 1.57 than at a Mach number of 1.96 in the frequency range of the flutter test 
results. The difference in overall levels is of the order that would be expected for this 
Mach number range in turbulent flow. (See ref. 20, for example.) 

The minimum Reynolds number, based on the distance from the leading edge of the 
splitter plate to the leading edge of the panel and minimum stagnation pressure, was 
approximately 1.3 X 10® at a Mach number of 1.96. The boundary -layer thickness ranged 
from 0.8 inch (2.03 cm) to 1.75 inches (4.445 cm) for the two Mach numbers at stagnation 
pressures of 718 psf (34.378 kN/m2) and 2160 psf (103.42 kN/m^). 
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MODE -SHAPE MEASURING APPARATUS AND PROCEDURES 

The vibration data measurements were made in the 55 -foot (16.764-m) vacuum 
chamber of the Langley dynamics research laboratory. In order to measure panel mode 
shapes remotely, a system was designed that incorporated a means of systematically 
handling a large volume of data. 


Apparatus 

A photograph of the splitter plate and the mode-shape measuring apparatus is 
shown in figure 18. Three 2-pound (8.90 N) shakers were positioned so that the first 
six streamwise bending vibration modes, with one half-wave in the cross-stream direc- 
tion, could be excited. A small length of compliant plastic tubing was attached to the 
shaker driving spindle. The shakers were mounted so that a 0.003 -inch (0.00762 cm) 
gap existed between the plastic tubing and the panel. When the panel was vibrated in a 
given mode, the spindles of the coils of the two inactive shakers were drawn back from 
the panel by a dc bias voltage on the shaker coils. The carriage shown in figure 18 
transported a displacement detector either continuously or incrementally by remote con- 
trol by use of two independent coordinate drive systems. The variable reluctance detec- 
tor was automatically maintained at 0.02 inch (0.0508 cm) from the panel surface. The 
ac analog signal from the detector was converted to a dc voltage for readout on a 
x-y plotter and a digital voltmeter. In order to minimize error caused by electronic 
noise, an error null circuit was used at the output of the ac-dc converter. This circuit 
provided a dc voltage equal in amplitude and opposite in polarity to the electronic noise 
signal. 

A reference displacement detector provided a continuous record of surface displace- 
ment at a point for the period of time necessary to map the mode shape of the upper half 
of the panel and thus provided the reference needed to correct the measured data for any 
variation in response amplitude with time. A reference detector was also used to pro- 
vide a response phase relationship over the vibrating surface. 

Procedure 

The vacuum chamber was evacuated to a density equal to the static air density of a 
point on the flutter boundary and the corresponding inplane load was applied. After a pre- 
liminary mode -shape search was made in the continuous mode of the carriage operation, 
data were taken in the incremental mode of operation in which the detector moved in pre- 
set increments over the panel. The amplitude, reference amplitude, phase, and surface 
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coordinates were automatically punched on cards at each point. The amplitude, phase, 
and one coordinate were also plotted simultaneously on the x-y plotter for visual 
reference. 

A program using Fourier series analysis was used to fit a curve through the data. 

A typical example is shown in figure 19 where the large symbols are the measured data 
and the small dots were obtained from the curve-fitting program. 

There were practical limitations to the mapping of the mode shapes of a panel under 
inplane load. The region of stress above the coalescence of the first two modes presents 
a problem because of the number of modes within a narrow range of frequencies of a 
damped system. Also, a very slight change in stress, due to small changes in tempera- 
ture for instance, caused the panel response to shift mode at a given frequency. As may 
be observed by referring to figure 7, the problem becomes increasingly severe with 
increasing length -width ratio. A system limitation was the inability to excite the first 
two modes at a length-width ratio of 3.77 and modes above the third for a length-width 
ratio of 1.0. It might have been possible to excite these modes if the system had been 
designed for a multiple shaker input in which a system in which the force is proportional 
to the deflection allows the use of multiple shakers without mode-shape distortion. 

As may be apparent from the foregoing discussion, the mapping of a sufficient num- 
ber of panel modes over a wide range of stress conditions with sufficient accuracy for 
use in flutter calculations presents a formidable task even for one panel and was found to 
be impractical. The system was useful for comparison of experimental mode shapes 
and frequencies with computed values and the determination of panel damping for the 
conditions of the flutter experiments. 
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DETAILS OF FLUTTER SOLUTION 
Zero-Stress-Mode Solution 

Since the flutter frequency appears as l/c^ in the off-diagonal as well as the 
diagonal elements in equation (7), the determinant of the coefficients was first put into a 
lambda form in which the eigenvalue, X = l/cu^, appears only in the diagonal elements. 

The determinant of the coefficients of equation (7) is in the form 

|a - cx| = 0 (Cl) 


where A and C are complex. The determinant is in the desired form when 


D - IX| = 0 


(C2) 


where 


[d] = 03 _1 Ca3 (cs) 

To invert the matrix [cf] by using real arithmetic procedures, the elements are sepa- 
rated into their real and imaginary parts. Rewriting equation (C3) in this form yields: 



(C4) 


The bar denotes a submatrix and here the Greek letters are used to denote the 
imaginary part of the complex matrix. Multiplying the terms and the left-hand side of 
equation (C4), and then equating the real and imaginary parts results in a matrix equation 
that yields the real and imaginary parts of 



(C5) 


The submatrices of equation (C5) expressed in the terms of equation (7) are: 
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In matrix form these equations become, for a four-mode solution: 
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(C7a) 


(C7b) 
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where 


D pg( m ) - Cmn 

(Lengthwise symmetric modes; 

m is odd) 

Dp A (m) = 

(Lengthwise asymmetric modes; 

m is even) 

c p SA ( m , n ) = c$ n 
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where 


X 'X 9 

yi l 


- N v C 1 

TOT - £ 

y 


c ss( m ) “ c mn 

(m = n, odd) 

(4) 

Css2( m > n ) = Cjnn 

(m / n, odd) 

C AA (m) = 

(m = n, even) 

^AA2^ m,n ^ = f 'mn 

(m ^ n, even) 


The mode shapes used in the calculations from equation (18) of reference 7 are: 


Xg = cos 2a £~ 


cos a 
cosh j3 


cosh 2/3£ (Symmetric modes) 




(C8) 


X A = sin 2 a£ - a - sinh 2/3£ (Asymmetric modes) 
A sinh 0 


The integers in equations (C7) designate the mode number of the arguments a 
and /3 which are defined in reference 7. Coefficients resulting from the integration of 
the products of the mode shapes over the length of the panel are given below. 


■>*<-> - w - 2 ^ <0, * , 


D PA (m) " S Rl(°m) ' 


2 D (°Tn>£m) 

C C (%i>^m) 
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C PAS (m ’ n) = ” A L1 ( a m ,Q h) + 
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A LLl(^^m) = °h (ftn shl °n cosh ^m ' cos % sinh &n) 
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The lambda determinant (eq. (C2)) was solved by means of a QR algorit] 
solution for a given panel, stress ratio, and air density was the value of k/^ 
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caused the imaginary frequency to be zero. The value of air density used in the calcula- 
tions corresponded to the experimental density. 

At low stress levels where there was no low-frequency flutter data, the initial value 
of dynamic pressure, and hence air density, was determined from an extrapolation of the 
experimental low-frequency flutter boundary. At these stress levels, the calculated 
flutter dynamic pressure was often higher than the dynamic pressure of the extrapola- 
tion. It was found that increasing the density in the calculations resulted in decreasing 
difference between the input and calculated density. Since the experimental data were 
obtained under variable density and constant velocity conditions, the appropriate flutter 
boundary was considered to be that which resulted in a constant flutter velocity. 


Stressed-Mode Solution 

The stressed-mode shapes and frequencies used in equation (12) were calculated by 
using the equations of reference 7. As was briefly discussed in the section dealing with 
panel frequencies, the equations fall into two regions, depending on panel geometry and 
inplane load. The parameter that governs which mode shape equations are used for a 
given mode of an isotropic panel is 


where 


( 0 > m \ 2 /n 2 /l_ £2 N x w2 cA 

\“r) ' W \y4C 0 + ^2 D C 0 ) 



(Cll) 


The value of B for a given mode decreases with increasing midplane compres- 
sive stress, and for the lower modes, becomes negative. The mode-shape equations 
(ref. 7, eqs. (18) and (19)) for B positive and negative (eqs. (C8)) are 


X Q = cos 2a ^ - 


cos a 
cosh /3 

sin a 


cosh 2/3£ (Symmetrical; B positive) 


X A = sin 2a £ sinh 2/3£ (Asymmetrical; B positive) 

A sinh /3 


(C12a) 


Xg = cos 2a ^ - — — — — ^ cos 2/3i~ (Symmetrical; B negative) 
cos /3 




X. = sin 2a ^ - sin a sin 2/3£ (Asymmetrical; B negative) 
A sin /3 


(C12b) 


32 



APPENDIX C 


The arguments, a and 0 or j3, and the frequencies were calculated for panels 1 
to 6 for stress-ratio increments of 0.10 between 0 and 0.30 and in increments of 0.05 
from 0.30 to 1.0. 

A sample determinant is given for which B is negative for the first two modes. 
The terms in this determinant would correspond to the first 4 by 4 terms of panel 2 at a 
stress ratio of 0.60. 


D r (D - ft 
K(2) C NAS (2 ^> 


K < 4 > CpNAS^ 4 ’ 1 ) 


DpA<4) 


K (D C nsa (1,2) 

0 

K(l) 


D ns (1) 

d nsM 


Dr (2) - £2 

K(2) C npas (2,3) 
d na< 2 > 


0 

= 0 (Cl 3) 

K(3) C pNgA (3,2) 
D ps( 3 ) 

Dr (3) - 

K (3) C pSA (3,4) 
Dp S (3) 


0 

K(4) C pA g(4,3) 
d pa( 4 ) 

D R (4) ~ ft 


D R (m) = 

( i ° r ) 2 (l 1 

*> 

\ 



\ m/ \ ^l k ^m 

) 



£ 

ii 

eR 

2 

1 (1 + ig) 


> 

(C14) 

K(m) = - 

( u r\ 2 1 




\“m/ 

J 



D NS< m > = C mn 

(m = 1,3,5,...; B 

negative) 

(Cl 5a) 

°NA^ m ) = c mn 

(m = 2,4,6,...; B 

negative) 

(Cl 5b) 

D p S ( m ) = c mn 

(m = 1,3,5,...; B 

positive) 

(Cl 5c) 

D PA M = eg), 

(m = 2,4,6,...; B 

positive) 

(C15d) 
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C NSA^ m,n ^ ~ C mn (m = 1,3,5,...; n = 2,4,6,...; B negative) 

(5) — 

Cnas^ 11 ) = c mn (m = 2,4,6,...; n= 1,3,5,...; B negative) 


C PSA( m,n ) = c mn (m = 1,3,5,...; n= 2,4,6,...; B positive) 


CpA S ( m > n ) = c mn (m = 2,4,6,...; n = 1,3,5,...; B positive) 


m = 1,3,5,...; B negative' 
positive j 


r (m — r< 5 ) ( m ~ B 

C NPSA (m,n) - °mn = g 

/c\ /m = 2,4,6,...; B 

^NPAS^ m,n ^ = I - 

\n = 1,3,5,...; B 


^PNSA^ m,n ^ = ^mn 


^PNAS^ m,n ^ _ ^rnri 


mn 


m - 2,4,6,...; B negative' 
positive; 

(m = 1,3,5,...; B positive ^ 
^ n = 2,4,6,...; B negative; 

= 2,4,6,...; B positive N 
^n = 1,3,5,...; B negative; 


(C15e) 

(C15f) 

(C15g) 

(C15h) 

(C15i) 

(C15j) 

(Cl 5k) 

(Cl 51) 


The code used to identify the integrals in equations (C13) and (C15) may be explained 
by considering an example. The integral ^pNSA^ m,n ^ * s P roc * uc t m th sym- 

metric (S) mode shape with B positive (P) times the first derivative with respect to £ 
of the nth asymmetric (A) mode shape with B negative (N). 

The integrals for equations (C15) are 

= ®R2 + ®Rs(^m) ^( a! m>^m) (C16a) 

^NA^ = ®Rl( CE ta) “ + H ^( Q! m>/ 3 m) (C16b) 

C NSA^ m,n ^ = A L2 ( Q5 ai ,CE n) " A L4( Q a , ^m) 

” B A4( CK ai , ^a) + F (°h»^n) ^(^l^m) B A2(^m>^n) 

(C16c) 
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^NAS^ m,n ) A Ll( a ni ,a n) + ^(^m^m) A L3( 0 ri , ^ m ) + ®A3( fIm, ^ n ) 


“ E (°h»0n) ^{ <x m*^m) B Al(^ m, ^ n ) 


^NPSA^ m,n ^ "" A L2( a ni ,a h) " G (°W0m) A L4( <3 h , ^ m ) 

2 c ( Q n>/ 3 n) B AA2 ) + 2 c ( Q b>/ 3 n) ^(^tn^m) B AA5(^m>^n) 
B BB( 0? m>^n) ■ D D(^m>/ 3 n) 


C NPAS^ m,n ^ A Ll( <x m» 0 ^) + H ( (X m»^m) A L3(°h ,/3 m) 


2 A A( Q h , £ } n) B AAl( 0em, ^ n ) + 2 A A( fl h , &0 H ( Q tn»^m) B AA6(^ni>^n 


B BB( 0?m, ^n) 


C PNSA^ m,n ^ ' A L2 (“m^) 


2 B B( 0! m>/ 3 m) A LL2 
A AA(°^ ,/3m ) 


- F ( Q h>/ 3 n) B A4 ( Q m , ^n) 


+ 2 F (°Vi>^n) B B( a m , ^m) B AA4 ( £*n ) 

B D(^m>^n) 


C PNAS^ m,n ^ ' “ A Ll(“m ,Q! n) 


2 D(a? m ,£ m ) A L.Ll( Q h , ^ m ) 


L AA 




+ E (°h>^n) B A3( (X m , ^ 3 n) 


2 E (%>^n) D (%i^m) B AA3 ( ^ m»^n) 

D D(^mA) 


The remaining integrals D pg , D pA , C p g A , C pAg may be found in equa 
tions (C9). The terms used in equations (C16) are: 




^m - 0m 


(C16d) 


(C16e) 


(C16f) 


(C16g) 


(C16h) 


(C17a) 
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^ L L3( QJ n ,/3 m) 

B A1 : 

B A2 i^m’ Pn) '' 

B A3( a?H1, ^ n ) : 
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, v cos 

E K’M%o Sftl 

(C17m) 

sin Ob 
■^(^h'^n) s ji n ^ 

(C17n) 

COS Obl 

G ( a m»^m) “ cos 

(C17o) 

sin /3m 

(C17p) 

= + ^n 

(C17q) 


The remaining terms Aj^, Aj^j ' ’ * ma y ^ oun( ^ equations (CIO). 

The determinant (eq. (C13)) was solved by the QR algorithm. In this solution, 
dynamic pressure and flutter frequency are obtained directly as a function of the struc- 
tural damping coefficient g. 
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TABLE I.- PANEL PROPERTIES 

[Panel, 2024-T3 aluminum; stiffener, 6061-T6 aluminum; frame, 6061-T6 aluminum] 



Length Width’ 


in. m 

13.5 

24.5 

24.5 
24.5 
24.5 

24.5 

13.5 

24.5 
24.5 
24.5 
24.5 
24.5 
24.5 

^Distance between rivet lines 


Panel 

Length thickness 



frame 


Panel 
flexural 
rigidity, D 


62 
122.488 
122.488 

62.714 
32.109 
32.109 

122.488 

256.75 

122.488 

62.714 

62.714 
62.714 
62.714 



Adjusted torsional stiffness, M 


Longitudinal 

stiffener 


lb-in^ 

0.8928 X 10 f 
1.2997 
1.2997 
1.1544 
.8529 
.8529 
1.0465 
1.4663 
.8048 
.7462 
.8928 
.8928 
.5229 


kN-m^ 

2.562 

3.730 

3.730 

3.313 

2.562 

2.448 

3.003 

4.208 

2.310 

2.141 

2.562 

2.562 

1.501 



Calculated 
(ref. 7) 


lb-in^ 


1.527 x lb 5 

1.735 

1.735 

1.527 

1.350 

1.350 

1.735 

1.959 

1.735 

1.527 

1.527 

1.527 

1.527 


psi kN/m^ 


Buckling stress, a cr 


Experiment 


Leading edge Trailing edge 


psi kN/m^ psi 


850 5 861 850 

1360 9 377 1240 

1770 12 204 1860 

1000 6 895 900 


jN/m* 

5 861 
8 550 
2 825 

6 206 


0 963 1000 6 895 900 6 206 

1702 11 735 1800 12 411 1800 12 411 

2101 14 486 2340 16 134 2160 14 893 

1315 9 067 1220 8 412 965 6 654 

2021 13 935 1930 13 307 1885 12 997 

2423 16 707 2460 16 962 2267 15 631 

2660 18 341 2780 19 168 2630 18 134 

1580 10 894 1050 7 240 1117 7 702 

2660 18 341 2220 15 307 2090 14 411 

1580 10 894 1275 8 791 1275 8 791 


























TABLE H.- EXPERIMENTAL AND CALCULATED TORSIONAL STIFFNESS 


Test 

specimen 

Design 

rivet 

spacing 

Thickness 

M t, 9 

Hat section 

Skin 

Calculated 

Percent difference 

= 100(1 - exp/calc) 

in. 

cm 

in. 

cm 

in. 

cm 

lb-in 2 

kN-m 2 

Design spacing 

1/2 rivet spacing 

1 

0.75 

1.905 

0.064 

0.163 

0.032 

0.081 

1.0922 x 10 5 

3.134 

22 

12.3 

2 

1.0 

2.54 

.090 

.229 

.050 

.127 

1.70 

4.878 

27.7 


3 

.75 

1.905 

.064 

.163 

.032 

.081 

1.0922 

3.134 

24 


4 

1.0 

2.54 

.050 

.127 

.040 

.102 

.975 

2.798 

21.1 

2.7 

5 

1.0 

2.54 

.032 

.081 

.040 

.102 

.716 

2.055 

20 

8.4 

6 

.75 

1.905 

.032 

.081 

.032 

.081 

.672 

1.928 

25 

11.7 

*7 

1.0 

2.54 

.090 

.229 

.040 

.102 

2.027 

5.817 

25 

15 


*Strap 0.10 in. (0.254 cm) on top of hat section to simulate panel frame. 


41 















TABLE HI.- VIBRATION TEST RESULTS 


Mode 

Structural 

damping 

coefficient 

Frequency, 

Hz 

°/<7 C r 

Pressure, 
mm Hg 

Panel 1 

1 

0.0296 

89 

0 

56 

1 

.0225 

67.5 

.60 

56 

1 

.0258 

90 

0 

143 

1 

.0233 

80 

.27 

143 

1 

.0474 

91 

0 

296 

2 

.0225 

150 

0 

56 

2 

.0283 

93 

.60 

56 

2 

.0295 

146 

0 

143 

2 

.0266 

116 

.27 

143 

2 

.0257 

146 

0 

296 

3 

.0292 

192 

.60 

56 

3 

.0338 

207 

.27 

143 

3 

.0328 

271 

0 

296 

Panel 2 

1 

0.0300 

69 

0 

63 

1 

.0320 

61 

.27 

63 

1 

.038 

75 

0 

125 

1 

.054 

82.4 

0 

305 

2 

.028 

92.5 

0 

25 

2 

.038 

56 

.45 

25 

2 

.0223 

93 

0 

63 

2 

.034 

78 

.27 

63 

2 

.0398 

93 

0 

125 

2 

.0300 

94 

0 

305 

3 

.020 

129 

0 

25 

3 

.036 

80 

.45 

25 

3 

.022 

120 

0 

63 

3 

.024 

104 

.27 

63 

3 

.022 

133 

0 

125 

3 

.022 

134 

0 

305 

4 

.020 

191 

0 

25 

4 

.0222 

130 

.27 

25 

4 

.032 

188 

0 

63 

4 

.036 

198 

0 

305 
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TABLE HI.- VIBRATION TEST RESULTS - Continued 


Mode 

Structural 

damping 

coefficient 

Frequency, 

Hz 

°/°cr 

Pressure, 
mm Hg 



Panel 2 



,5 

0.008 

263 

0 

25 

5 

.036 

200 

.45 

25 

5 

.020 

233 

.27 

63 

5 

.010 

262 

0 

125 

5 

.012 

260 

0 

305 



Panel 3 



1 

0.0121 

94 

0 

1.0 

1 

.032 

98 

0 

105 

1 

.032 

92 

.45 

105 

1 

.0618 

107 

0 

185 

1 

.044 

100 

.20 

185 

1 

.056 

104 

0 

250 

2 

.010 

131 

0 

1.0 

2 

.02 

131 

0 

105 

2 

.02 

106 

.45 

105 

2 

.02 

125 

0 

185 

2 

.02 

118 

.20 

185 

2 

.022 

128 

0 

250 

3 

.018 

176 

0 

1.0 

3 

.020 

176 

0 

105 

3 

.0385 

165 

0 

250 

4 

.026 

225 

0 

1.0 

4 

.030 

225 

0 

105 

4 

.0362 

166 

.45 

105 

4 

.0280 

203 

.20 

185 

4 

.025 

222 

0 

250 

5 

.018 

310 

0 

1.0 

5 

.040 

314 

0 

105 

5 

.024 

272 

.45 

185 

5 

.0142 

296 

0 

250 

Panel 4 

1 

0.0221 

105 

0 

44 

1 

.0380 

104 

.55 

44 

1 

.0227 

107.5 

0 

116 

1 

.0210 

108 

.20 

116 

1 

.0467 

110 

0 

286 


TABLE m.- VIBRATION TEST RESULTS - Concluded 


Mode 

Structural 

damping 

coefficient 

Frequency, 

Hz 

%cr 

Pressure, 
mm Hg 

Panel 4 

2 

0.033 

135 

0 

44 

2 

.0411 

128 

.55 

44 

2 

.0345 

133 

0 

116 

2 

.0309 

132 

.20 

116 

2 

.0322 

128 

0 

286 

3 

.0166 

153 

0 

44 

3 

.0332 

123 

.55 

44 

3 

.0190 

150 

0 

116 

3 

.0236 

140 

.20 

116 

3 

.0160 

150 

0 

286 

4 

.0168 

194 

0 

44 

4 

.0196 

138 

.55 

44 

4 

.0203 

192 

0 

116 

4 

.0196 

172 

.20 

116 

4 

.0320 

194 

0 

286 

5 

.0223 

243 

0 

44 

5 

.0495 

243 

0 

286 

Panel 5 

2 

0.0456 

183 

0.45 

103 

3 

.0260 

195 

0 

9 

3 

.064 

190 

0 

103 

3 

.040 

181 

.45 

103 

3 

.062 

192 

0 

200 

3 

.0358 

186 

.27 

200 

3 

.0300 

192 

0 

300 

4 

.03 

227 

0 

9 

4 

.056 

227 

0 

103 

4 

.0600 

242 

0 

200 

4 

.024 

200 

.27 

200 

4 

.03 

232 

0 

300 

5 

.020 

264 

0 

9 

5 

.036 

269 

0 

103 

5 

.042 

280 

0 

200 

5 

.020 

228 

.27 

200 

6 

.02 

314 

0 

9 

6 

.024 

225 

.45 

103 


44 




Figure 1.- Panel mounted on splitter plate. L-65- 1693-1 

cn 




Figure 2.- Splitter plate cavity. 


L-65-1876 .1 







Figure 3.- Under side of panel 6. 


L-65- 7638.1 


-a 




^ 0 . 90 - 

( 2 . 29 ) 






L. E. 


a, b, c 





(a) Panel 1; Z/w = 1.0; (f) Panel 6; Z/w » 4.202; (J) Panel 10; Z/w - 3.77; 

h = 0.040(0.1016); h = 0.032(0.0813); h = 0.040(0.1016); 

P P P 

h = 0.064(0.1626) h = 0.070(0.1778) h = 0.050(0.127) 

s s s 


h = 0.040(0.1016); 
P 


h = 0.050 (0.127) 

s 





(b) Panel 2; Z/w = 1.96; (g) Panel 7; Z/w = 1.0; (k) Panel 1 1 ; Z/w = 2.88; 

h = 0.050(0.127); h = 0.050 (0.127); h = 0.040 (0.1016); 

P P P 


h s =0.080(0.2032) 


h s = 0.070(0. 1778) 


h = 0.064(0.1626) 
s 





(c) Panel 3; Z/w ■ 2.43; (h) Panel 8; Z/w = 1.96; (I) Panel 12; Z/w = 3.77; 

h = 0.050(0.127); h =0.064(0.1626); h =0.040(0.1016); 

P P P 

h s = 0.080(0.2032) h s = 0.080(0.2032) h s = 0.064(0.1626) 


h s = 0.064(0.1626) 





(d) Panel 4; Z/w = 2.88; (i) Panel 9; Z/w = 2.88; (m ) Panel 13; Z/w ■ 2.88; 
h = 0.040(0.1016); Z/w = 2.88; 

P » / » h = 0.040(0.1016); 

Q ' 080 (Q-2032 ) h p = 0.050(0.127); P __ Q ^ {Q Qm) 

KZZZZZ2 h s = 0.050(0.127) 


(e) Panel 5; Z/w = 3. 7T; 

h = 0.032(0.0813); 
P * 

b =0.080(0.2032) 


on panel denotes calculated natural frequencies, 
dynamic pressure at flutter and flutter 
frequencies 

on panel denotes measured natural frequencies 
and structural damping 

on panel denotes panels adjacent to center test 
bay are 0.040 (0.1016) 


Figure 5.- Summary of panel configurations used in investigation. Thickness of flutter panel h p and primary longitudinal stiffener h s 

given in inches (cm). 
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(a) Longitudinal stiffeners. (b) Transverse frame. 

Figure 6.- Examples of torsional stiffness test specimens. 
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Calculated 

hyperbolic-trigonometric solution 


Trigonometric solution 


Experiment 


0.4 0.6 

°"/ r 


(b) Panel 2 ; l/w = 1.9t 
Figure 7.- Continued. 
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Frequency, Hz 








f - 108 Hz — - 


f = 113.2 Hz 


(a) cr/cr,^. “ < 
I mm Hg 


f *101 Hz "A . - 


f = 110.2 Hz 



f = 148 Hz 



(b) °* 20 °‘> 
187 mm Hg 


f = 93 Hz 


Calculated 

Measured 


f = 106.1 Hz 





(c) cr/cr. r = 0.450; 
105 mm Hg 


Figure 8.- First- and third-mode shapes of panel 3 after flutter tests. 
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Minimum flutter dynamic pressure, 
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Figure 12,- 





Flutter dynamic pressure, q, psf Flutter frequency, Hz 





Flutter dynamic pressure, 
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Id) Panel A 


Figure 14. 


o 



0.6 0.8 1.0 1.2 

cr /°c r 


i; l/w = 2.88. 


Continued. 




Flutter dynamic pressure, q, psf Flutter frequency, Hz 


I 




(e) Panel 5; l/w = 3.77. 
Figure 14.- Continued. 
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I 


Flutter dynamic pressure, 


Flutter dynamic pressure, q, psf Flutter frequency, Hz 



3000 


2500 H 


2000 


1500 h 


1000 h 


500 h 


Calculated 

Zero-stress-mode shapes 
Zero-stress-mode shapes 
Stressed-mode shapes 
O Experiment 
# No flutter 
© Deeper penetration 



H 120 


Hi oo 


H 80 


H 60 


H 20 


cr 



(f) Panel 6; l/w = 4.202. 
Figure 14.- Concluded. 
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Figure 17.- Pressure spectra and overall fluctuating pressure normalized by dynamic pressure at Mi = 1.57 and M 



Figure 18.- Mode shape survey apparatus in Langley 55-foot vacuum chamber, 


L-66-339.1 


8 


Measured value 


Fourier series interpolation 


• • 



I 


-2k- 


-4h 


-6 k 


“8h 


-10 


0 




.4 


« 







Jl 

.6 



J 

1.0 


Figure 19.- Typical measured mode shape and program curve fit. 
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